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Abstract 

We study Zjv strings in nonabelian gauge theories, when they can be considered 
as domain walls compactified on a cylinder and stabilized by the flux inside. To make 
the wall vortex approximation reliable, we must take the 't Hooft large N limit. Our 
construction has many points in common with the phenomenological bag models of 
hadrons. 
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1 Introduction 

In a previous paper P], we explored the idea that a flux tube can be though, 
under particular circumstances, as a domain wall compactified on a cylinder. To 
make it possible, the theory must admit two degenerate vacua: one in the Coulomb 
phase, that is kept inside the cylinder, and the other one in the Higgs (or confining) 
phase, that is kept outside the cylinder. For the simplest example, the abelian Higgs- 
Coulomb model, we were able to determine a condition under which the wall vortex 
approximation is quantitative reliable. When we increase the number n of quanta of 

2 

magnetic flux, the radius of the vortex grows like ?7,3, while the thickness of the wall 
remains fixed. When n is enough large, the thickness of the wall becomes negligible 
and the tension of the vortex is given by the simple minimization of the energy density. 

The present paper started with the following question: can we find a nonabelian 
version of the wall vortex? In P we found examples where Z^v strings in a confining 
SU{N) gauge theory can be qualitative thought as wall vortices. What we are locking 
for, is some realization in which the wall vortex in quantitative reliable, so that we 
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can apply the simple argument of the energy minimization used for the abelian Higgs- 
Coulomb model. 

If we want to mimic what we obtained for the abelian model, we must increase the 
number n of flux quanta, until the radius of the n-string becomes much larger than 
the thickness of the wall. A problem immediately arise: in a SU{N) gauge theory, 
n is limited to be smaller than N. This means that, to have a chance of obtaining 
some result, we must explore the large limit of the nonabelian gauge theory j2]. 

In the first part of the paper we will explore the case of degenerate vacua, one in 
the Coulomb phase and the other one in the confining phase, and we will consider 
two examples of supersymmetric gauge theories. In the second part of the paper we 
face the more general situation in which the Coulomb vacuum is metastable (or in 
the extreme case instable). Our hope is to apply these ideas to non supersymmetric 
gauge theories and we will try to do that for large N QCD. 

The ideas we will expose are not new to particle physics. Almost thirty years ago 
a lot of works have been made on the bag models of hadrons. The first one was the so 
called MIT bag model [SI, where the quarks where modelled as free fermions bounded 
in region of finite volume and non zero energy density. This model had a great success 
in explaining static properties of hadrons j3]. A lot of different version of the bag 
model have been proposed after that. For example, in the so called SLAC bag model 
jS], the bag had finite tension and the interior of the bag had zero energy density. 
Another interesting approach was the Friedberg-Lee model of hadrons In this 
works, by means of an auxiliary scalar field, they derived the bag as a domain wall 
interpolating between a metastable vacuum and a true vacuum. In the metastable 
vacuum the quarks have a small mass, while in the true vacuum they have a great (or 
infinite) mass. In the Friedberg-Lee model, hadrons where realized as nontopological 
solitons. Despite their success in explaining static properties of hadrons, bag model 
had some difficulties in describing interactions and so they where though only as 
phenomenological modelling of the real theory, QCD. Attempts to derive the bag 
model from QCD first principles have been made in j3| and [Hj. 

In this paper we will argue that some confining gauge theory can share, in the 
large N limit, some properties in common with bag models. To be more clear, we give 
now a definition of what we mean by a bag model. The definition is not precise but 
sufficiently large to include all the possible realizations. In a bag model we have two 
phases: one in the interior of the bag that can contain quarks and gluons, the other 
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in the exterior of the bag where only gauge singlets can live. The interior Coulomb 
phase can have some energy density Eq and the bag can have a tension Tw- For 
example, in the MIT bag model Sq 7^ and Tw = 0, while in the SLAC bag model 
£o = and Tw 7^ 0. In general we can define a bag model action: 

L = [ d?xC,nt-eoV -Tw I d^e V- det h , (1.1) 

where Cint describes the dynamics in the interior Coulomb phase. Another important 
ingredient, to complete the definition of the theory, consist in specifying the boundary 
conditions of the fields at the bag surface. This condition is that color must be trapped 
inside the bag. In Figure Q there is a meson in the bag model. When the meson is 
rotating very fast, it becomes approximatively the wall vortex with a quark and an 
antiquark at the ends. The connection between the bag and the string was first 
recognized in ^j. 

Confinig Phase 




Figure 1: A meson in the bag model. Quarks and gauge fields are in a Coulomb phase trapped 
inside the bag. The Coulomb phase has energy density Eq and the bag has surface tension Tw- The 
confining vacuum is outside the bag. 

In the large limit of a SU{N) gauge theory, Feymann graphs organize them- 
selves into a genus expansion in powers if i 0. For this reason it is believed that 
some dual string theory should describe nonabelian gauge theories, and this string 
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theory should become weakly coupled as goes to infinity. The AdS/CFT corre- 
spondence is a concrete realization of this duality ^Oj for the A/" = 4 superconformal 
gauge theory. After that, some examples have been found where the gauge theory is 
not superconformal but confining fJEj- For the above mentioned reasons it may 
be sound a bit strange that a confining gauge theory becomes well approximated by 
a bag model in the large A^ limit. An important point, is that our claims regard 
n-strings with n of order A^, and not the fundamental string over which is supposed 
to be built the dual string theory. So the hadrons that becomes a bag are the exotic 
mesons g"-^" with n N } 

The string tensions Ty{n) has been much investigated, both theoretically and 
numerically (see fl] for a review), in particular the ratio of string tensions defined 
by A^) = . From the theoretical side there are two important predictions: 
the Casimir scaling and the sine formula: 

Casimir scaling. In an intermediate range of distance, the force between two static 
sources is proportional to the quadratic Casimir of the representation [T5t IT6]. 
As A^ becomes large, this range goes to infinity and the extrapolation suggest 
the ratio for string tensions 

niN — n) , , 

Rin,N)= (1.2) 

This procedure is not rigorous, since we should first make the distance to go to 
infinity, and then make the large A^ limit. It has be noted that the Casimir 
scaling, when n is kept fixed and A^ goes to infinity, has corrections jj instead 
of the expected This should rule out the Casimir scaling but there are also 
different opinions on that 



Sine formula. This formula first appeared first in [20] 

R{n,N) = '^^, (1.3) 
sm (^) 

studying M = 2 SU (N) gauge theory, softly broken by the adjoint mass term 
fiTi $2. In |2I!, in the MQCD contest, has been shown that this formula is true 
also in the opposite limit fi ^ oo. In j22| it has been shown that the sine formula 



^This should also overcome previous problems founded in trying to relate large N QCD and bag 
models ^3). 
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doesn't hold in the intermediate regime but has non universal corrections. The 
formula reappeared in j2Hj in the contest of cascading gauge theories. Another 
meaning of the sine formula has appeared recently in [21]. The conclusion is 
that the sine formula, even if not directly derived in the QCD contest, present 
some universal characteristics that makes it a good candidate for QCD, or for 
any confining gauge theory. 

On the lattice side, a lot of works have been done to compute the tensions in pure 
SU{N) Yang-Mills ISniEHj- The most recent work on the subject gives results 
up to = 8. In the end of the paper we will confront our theory with these results. 

The paper is organized as follows. In Section |21 we consider the wall vortex in the 
abelian Higgs model. First we review the results of P in the case where the Coulomb 
phase has the same energy of the Higgs phase. Then we extend the results to the 
most general case where the Coulomb phase is a stationary point of the potential 
(this must always be due to the ^(1) symmetry). In Section El we study a model in 
which there are Z^r solitonic magnetic strings that in the large limit become wall 
vortices. In Section |3] we try to apply our ideas to confining string in ordinary 
QCD. Even if our reasoning doesn't follows directly from QCD first principles, it gives 
a sharp prediction for the ratio of string tensions in the large A^ limit. The future 
lattice computations should easily prove or disprove our hypothesis. 



2 The Abelian Theory and the Wall Vortex 

We consider the abelian gauge theory coupled to a charged scalar field 

^ = -^2F,uF^'' - m - ^A,)q\' - Vi\q\) . (2.1) 

In P we considered the case in which the potential has two degenerate vacua, one 
in the Coulomb phase and the other in the Higgs phase. In 12.11 we briefiy review the 
results that we have obtained. In 12.21 we make a step forward trying to generalize 
the wall vortex idea to the most general case, where the Coulomb phase is not a true 
vacuum but only a stationary point. 



5 



2.1 Coulomb-Higgs model and the surface bag 



Now we consider the potential of Figure |21 There are two degenerate vacua: one in 
the Higgs phase where \q\ = qq, and the other in the Coulomb phase. For this reason 
we call it abelian Coulomb-Higgs model. The magnetic vortex [SHI IHO] in the Higgs 



Vi\q\) 




qo \q 



Figure 2: A potential with two degenerate vacua: g = is in the Coulomb phase while \q\ ~ qo is 
in the Higgs phase. 

vacuum is nothing but the wall interpolating between the two vacua, compactified on 
a cylinder and with the Coulomb vacuum left inside. The energy density as function 
of the radius is 

^(^^ = 2^^ + ^^2^^ ' ^^-^^ 
where $b is the magnetic flux. The magnetic flux is quantized in integer values 
$B = 2nn. The stable configuration is the one that minimizes the tension: 



Ty = S^VT {^] \ Ry = r2{l^. (2.3) 



w 



In this simple calculation we have neglected the thickness of the wall A^i/ and this 
is the crucial point. We can trust ()2.3j) only if the radius if the vortex Ry is much 
greater than Avi/- In P we argued that ()2.3j) can be used in a self-consistent way to 

2 

determine its region of validity. Taking Ry = \f2 ^"^ i for true, we increasing n 
keeping fixed all the other parameters of the theory. Note that Ayj/ does not depend 
on n. There will be some value n* above which Ry ^ A^y and so ()2.3|) can be trusted. 
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2.2 Higgs model and the volume bag 

Now we discuss an abelian gauge theory with a potential hke Figure El The Coulomb 
vacuum is metastable and has energy density Eq- 




qo \q 



Figure 3: The vacuum g = is in the Coulomb phase, it is metastable and its energy density is 
Eq. The other vacuum at \q\ = go is in the Higgs phase and has zero energy density. When Eq = 
the potential becomes the one of Figure [3 and when £o = vq it becomes that of Figure 01 

Also in this case is convenient to think of the flux tube as a domain wall compact- 
ified in a cylinder, with the metastable Coulomb vacuum inside and the true Higgs 
vacuum outside. Neglecting the thickness of the wall, we can write the tension as 
function of the radius:^ 

T{R) = ^ + Tw27rR + EotiR^ . (2.4) 
There are two regimes in which ()2.4j] can be easily solved. 

Surface (or SLAC) bag. This region is when n satisfies the conditions 

< n < —— . (2.5) 



^To derive this formula we could first think of the vortex tension as function of two parameters, 
the radius R and the thickness A of the shell where the q field goes from to go- The tension is 
roughly T{R, A) = + + RAvq + e^R^ . Minimizing with respect to A we obtain A^/ ~ ^ . 
So we can interpret it as a domain wall since is independent on n and R. 
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In this limit the surface term in ()2.4|) dominates over the volume term and the 
minimization gives: 



Ty = S^TT (^] \ Rv = ^^H^. (2.6) 



w 



Note that the surface region ()2.5p exists only if Eq vq. As £o is increased until 
it reaches Vq, the SLAG region is eaten by the MIT region. 

Volume (or MIT) bag. This region is when n satisfies the condition 

2 

< n . (2.7) 



In this limit the volume term in ()2.4j) dominates over the surface term and the 
minimization gives: 



T^ = 2V27r^, Ry = 4^ U^. (2. 
Note that the tension is proportional to n, as happens in the BPS case. 



General conjecture 

The general conjecture is that the result (j2.8p works for every potential, even if 
the Coulomb phase is not metastable but instable like in Figure |3] Let's write the 
conjecture for clarity. Consider the abelian Higgs model h2.1\) with a general potential 
that has a true vacuum at \q\ = go 7^ and a Coulomb phase with energy density 
V{Q) = eo 7^ 0. Call Ty{n) the tension of the vortex with n units of magnetic flux. 
The claim is that 

lim Tv{n) = 2v^7r . (2.9) 

n— »oo e 

We will give a substantial check of this statement in |2H] using numerical computa- 
tions. 

A check 

Now we make a non trivial check of the conjecture using the famous example solved 
by Bogomoln'y jl^. When the potential is 

V{\q\)='^{\<P\'-0\ (2.10) 
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Figure 4: In the extreme case Eq = vq the Coulomb vacuum is still stationary but instable, 
the tension is 

Tv = 27in^ (2.11) 

for every value of n. Solving the model with our trick, the result must coincide with 
()2.11|) . For the BPS potential ()2.10p . the energy density of the Coulomb vacuum is 
£o = ^-2~ cind, using ()2.8|) . we find exactly ()2.11|) . This can hardly be considered only 
a coincidence. 



3 Solitonic Z^v Strings 

The principal aim of this paper is to see if it possible to have a nonabelian gener- 
alization of the wall vortex. In this section we consider a toy model of solitonic and 
magnetic strings. To have this we need an SU{N) gauge theory with scalar fields 
that are not charged under the center of the group. If in some vacuum the scalar vevs 
brake completely SU{N), then the solitonic strings are stabilized by the homotopic 
group 

In the following we will consider the a broken version of the A/" = 4 supersymmetric 
gauge theory. There are three adjoint scalar fields that break completely the gauge 
group and are responsible for the formation of the Ztv strings. 
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3.1 J\f =1* and the surface bag 



The first case that we consider is the A/" = 1* SU{N) gauge theory, that is A/" = 4 
broken to A/" = 1 by mass terms for the chiral superfields $1, $2 and $3. The 
superpotential is the A/" = 4 interaction plus the mass terms: 



WV=i* = ^Tr ($1 [$2, $3] + ^$1^ + ^$2^ + ^$3^) . (3.2) 



The stationary equations are, a part from numerical factors, the commutation rela- 
tions of the SU{2) algebra [HSEB]- For example, deriving ()3.2|1 with respect to $3, 
we obtain $2] = ~''^3$3- Making the rescalings 

$1 = iy/m^m^^i , $2 = iy/m^ml^2 , ^3 = i^/m^m^^^ , (3.3) 
we obtain exactly the SU{2) algebra as equations of motion: 



$. <T). 



ieijk^k ■ (3.4) 



The vacua of the theory are obtained choosing a partition of 

N 



J2dkd = N, (3.5) 



d=l 



SO that the N x N matrices $4 are covered with spin representations of the SU{2) 
algebra. The gauge group is classically broken 

SU{N) ^MM , (3.6) 
and the vacua of the theory are divided into two classes: 



Massive Vacua. These vacua are called massive because there is a mass gap. For 
every divisor of we must cover the matrices $i only with representation of 
the same dimensionality. In this case the gauge group is classically broken to 
SU ) and there are no f/(l) factors. The SU (f ) group confines and there 
is a mass gap. This vacua are in one to one correspondence with the possible 
phases of a general confining gauge theory . The works [321 El] showed that 
the SL{2, Z) duality of the original J\f = 4: theory, exchange the massive vacua 
among them. 
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Coulomb Vacua. When there are at least two representations with different dimen- 
sionality, the unbroken gauge group has at least one U{1) factor. Since there 
is no strong dynamics for the [/(l)'s, they survive in the infrared giving some 
massless particles. 

Before proceeding to the central point, we need to show some qualitative properties 
of the domain walls in the A/" = 1* theory. What is really important for us in the 

dependence of the tensions and of the thicknesses of the walls. If we assume 
that the walls are BPS saturated, their tension is governed by the difference of the 
superpotential in the two vacua Fo^^ ^ general vacuum ()3.5j) the superpotential 
()3.2j) is proportional to the sum of the Casimirs of the spin representations 

W r^m'^N^ (^^kd{d-l){d+l)^ , (3.7) 

where for simplicity we have used rh = ^fm^rrh^m^. And so the superpotential, goes 
like ~ A^^. Choosing two generic vacua, the interpolating wall scales like 

TvK-A^^m^, Av^~i. (3.8) 

m 

There are some exceptions to this scalings. Take for example the domain wall between 
the Higgs vacuum and the Coulomb vacuum where the is partitioned into a — 1 
and a 1 representation (call it (A^ — 1,1) for simplicity). The leading terms in the 
superpotential cancel each other and the wall goes like:^ 

Tw^N''m\ A^^^. (3.9) 

A* m 

Now we make our claim. The solitonic Z^v strings in the Higgs vacuum are made 
by a domain wall compactified on a cylinder with a Coulomb vacuum inside. When 
A^ becomes large the energetically favorite Coulomb vacuum is the one where A^ is 
partitioned into a A^ — 1 and a 1 spin representation. For sufficiently large A^ the 
radius of the n-strings becomes much larger that the thickness of the wall, and so 
that the wall vortex condition (namely ^ Rv) is satisfied. 

First of all we check that the wall vortex condition is satisfied. The Coulomb 

vacuum of interest, has only one U{1) factor, the one generated by the matrix 

■^This mechanism is very similar to what operates in A/" = 1 SYM where the tension of the 
domain waU between adjacent vacua goes hke instead of the expected |88| . This point has 
been clarified in and PU] , 
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diag (1, . . . , 1, — (A^ — 1)). This is the only generator in the Cartan subalgebra that, 
when exponentiated, passes through all the elements of the center of SU{N). In 
particular the n-string has charge 



n 
N 



( 1 



v 



\ 



-(iv-1) ; 



The tension as function of the radius R is 



and, minimizing with respect to i?, we obtain: 



\N ) m 



(3.10) 



(3.11) 



(3.12) 



Since A\y ^ jf, ioi N sufficiently large Ry is much greater than A^y and the wall 
vortex approximation works. The ratio of string tensions is thus: 



R{n, N) = min (na , (A^ — n) 3 ) . 



(3.13) 



Now we consider the other point: why the Coulomb vacuum (A^ — 1,1) should be 
preferred with respect to the others? Or again: why the strings should all choose 
the same Coulomb vacuum? We don't have a rigorous proof for these questions but 
only an argument in favor of it. For example consider the n-string. Its flux could also 
be generated by the U{1) of the Coulomb vacuum (A^ — n,n): 




-{N-n)lr 



(3.14) 



Note that this ^7(1) cannot reach all the elements of the center of SU{N) so, if this 
generator comes out to be energetically favorite, we would have a different U{1) inside 
every n-string. What we are going to do is to evaluate the tension and compare it 
with the one obtained with the Coulomb vacuum (A^ — 1,1). The tension as function 
of the radius is 



TiR) 



n(N — n) 



m^N^R . 



(3.15) 
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The flux term has an ""^ from the trace of the square of ()3.14|) . and a power 
from the 't Hooft scahng. The domain wall is an ordinary soliton that scales like 
()3.8j) . The minimization of ()3.15j) with respect to R gives: 

Ty^^n{N-n)N^m' , i?^, ^ . (3.16) 

V JS^ m 

2 

Note that in this case the radius scales like A^~3 if we send n to infinity like A^. So we 
can not apply our approximation since the thickness of the wall ()3.8j) remains finite. 
Note that (j3.16|) would give a tension that scales like would be greater than the 
other one ()3.12j) . But, as we said, the wall vortex approximation doesn't work and we 
have to find another way to estimate the tension. The minimization of ()3.16p would 

2 

give Rv ~ N~3 that is much lesser than the thickness of the wall Aw ~ 1- This 
means that the scalar fields are spread all over the radius of the vortex and so their 
kinetic energy goes like ^ and dominates over the magnetic energy term that goes 
like This imply that the vortex tension is essentially given by the minimization 
of the scalar field action A^^ {dcpdcp + V{(fi)) and so it scales like A^^. This shows that 
the (A^ — 1, 1) vacuum gives a lighter tension than the (A^ — n, n) vacuum in the large 
A^ limit. 



4 Large iV QCD 

Now we try to see if there is any chance that the wall vortex scenario is realized 
in large A^ QCD. 

First we recall some well established results regarding the QCD string tension 
at large A^. Consider the interaction between two fundamental strings. It has been 
shown in jTTj that interaction vanishes like Thus the tension of the n-string, when 
n is fixed and A^ becomes large, is equal to n times the tension of the fundamental 
string 

Ty{n) = nTy{l) + O 1^-^^ . (4.1) 

If instead we keep n of order A^, the interactions are of order 1. In fact, there are (2) 
ways of making a fundamental interaction between any couple of strings, and so in 
total we have an interaction of order i^^- 

Another information we will use is that the tension of the fundamental string 
remains of order 1 as A^ is increased. In fact its tension determines the mass of the 



13 



meson and, for large A^, these masses remains finite 

Now we explore the possibility that the Z^v strings of QCD becomes wall vortices 
in the large limit with kept fixed. A sort of magnetic dual of what happened in 
the = 1* theory, but with a Coulomb energy density different from zero. We write 
the tension as function of the radius of the vortex considering only the dependence 
on n and A^: 

T{R) -J^ + N-K'R' . (4.2) 

Note that the energy of the electric flux term is instead of the obtained for 
the solitonic magnetic vortex in ()3.1H) . The reason is this: normalizing the gauge 
kinetic term as we have done, ^FF, the magnetic monopoles have charge of order 

2 

1 while the electric particles have charge of order ^ and this brings the A^ factor in 
the denominator of the electric field energy density. The Coulomb energy density is 
instead N°'A^R'^ where is just for dimensional reasons and a is a parameter that 
we are going to fix. 

Minimizing ()4.2j) with respect to R we obtain 



n . o „ aI m? 1 



Tv ~ A^ , Rv ~ {/iTTT— T • (4-3) 

We could have a spectrum like ()4.H1 . only if the a parameter in ()4.2|1 would be equal 
to 1. This in fact is the only way to match ()4.3|) with ()4.1|) since 7V(1) is of order 
1. In this way the tension is of order nN^ and so, keeping fixed n and sending A^ to 
infinity, we obtain a finite limit. 

Now there are two left points to explain. Let us rewrite the energy density, the 
tension and the radius now that we have fixed a = 1: 

Fn 1 

eo~A^^A, Ty~r^A^ Ry ^ J - - . (4.4) 

First we have to explain how is possible to have a Coulomb energy density that 
goes like A^A^. Second we have to explain how the wall vortex approximation can 
work, that is we have to find a domain wall between the confining vacuum and the 
Coulomb vacuum with a thickness that goes like Aw ~ j^. Both of these points can 
be explained by an effective Lagrangian that scales with an overall factor of A^^ 

Ceff = N'F[B,dB,...] , (4.5) 
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and a distance between the confining vacuum and the Coulomb phase that is 6B ~ 
Essentially is the same mechanism that worked in Section IHl^ 

4.1 Lattice data and saturation 

Now we confront our theory with the lattice data. First of all we make a brief 
discussion of the saturation limit. This is the best way to confront the experimental 
results for different gauge groups. Given any ratio of string tensions R{n, N) = , 
we rewrite it as a function of the ratio x = j^, and then rescale with a factor. For 
example the Sine formula and the Casimir formula become respectively: 



Note that in this way we can plot j^R{n, N) with respect to x in the same graph for 
all the values of N. The saturation limit is when the corrections can be neglected. 
This also the limit in which the various data becomes dense and describe a continuous 
curve. The Sine formula saturates to ^ sin (nx) while the Casimir formula saturate to 
x{l — x). An important thing to note is that both the Sine and the Casimir formula 
saturate from above, that is the o{j^) corrections are positive . Our formula is instead 



Since the deviation from an exact wall vortex is given by the wall thickness of order 
is natural to have jj: correction. It's also natural to expect a negative correction 
and so a saturation from below. 

Finally we confront with the lattice experiments that are plotted Figure El Up 
to now the largest for which computations have been done is = 8. In principle 
these data could be consistent with the formula min (x, 1 — x). To explain the ob- 
served deviation we should have corrections with coefficients of order 1. When the 
saturation will be reached it will be easy to prove or disprove our formula. 

"^It's also the same mechanism that works for the pure J\f = 1 theory for domain walls between 
two adjacent vacua |H81 1391 1^ . The effective Lagrangian can be written as £e// = N^[^ + A]. If 
6B ^ jj, minimizing with respect to A we obtain A\y ^ i and Tw ^ 




(4.6) 




(4.7) 




(4.8) 
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Figure 5: In the three graph are reported the lattice data taken from [^El- They refer respectively 
to SU{4:), SU{6) and SU{8). The lines plotted are respectively the Casimir formula (green line), 
the Sine formula (red line) and the min(a;, 1 — x) formula (blue line). 
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